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ABSTRACT  :  -To  determine  the  noncentral  gravitational  field  of  the  Moon,  the 
motion  of  the  artificial  Moon  satellite  "Luna  10"  is  analyzed  with  gravitational 
forces  of  the  Earth  and  the  Sun  taken  into  account.  From  the  statistically 
processed  measurements  of  the  trajectory  of  Luna-10  carried  out  over  the  period 
of  its  existence  (from  April  3  to  May  30,  1966),  it  is  established  that  the 
effect  of  noncentrality  of  the  gravitational  field  of  the  Moon  is  an  essential 
factor  in  the  evolution  of  the  orbit  of  Luna-10.  The  perturbations  of  its 
circumlunar  orbit  due  to  noncentrality  of  the  gravitational  field  of  the  Moon 
are  particularly  noticeable  in  the  evolution  of  the  longitude  of  the  ascending 
node  St  and  the  angular  distance  of  the  pericenter  to  the  ascending  node 
The  perturbations  of  elements  -O^and  of  Luna-10  during  the  time  of  its  existence 
(460  revolutions)  due  to  noncentrality  of  the  gravitational  field  of  the  Moon 
and  due  to  gravitational  effects  of  the  Earth  and  the  Sun  (the  gravitational 
effects  of  other  planets  are  neglected)  are  presented  as  functions  of  the 
•number  of  revolutions  of  the  Moon  satellite.  A  comparison  of  the  results  shows 
that  perturbations  for  Luna-10  due  to  noncentrality  of  the  gravitational  field 
of  the  Moon  exceed  the  perturbations  due  to  gravitational  forces  of  the  Earth 
•and  the  Sun  5  to  6  times. (  >The  gravitational  potential  of  the  Moon  is  taken  in 
the  form  of  expansion  in  ^herical  functions  with  unknown  coefficients  (the 
gravitational  field  parameters)  and  d^  (n  =  2,  3,  .  .  .,  m  =  0,  1,  .  .  .n)  . 

As  a  result  of  processed  measurement  data,  the  numerical  values  for  eleven 
expansion  coefficients  are  derived  which  are  presented  with  maximum  possible 
errors.  To  illustrate  the  determined  gravitational  potential  of  the  Moon,  the 
level  surface  passing  through  the  point  with  spherical  coordinates  r  =  1738  km, 

=■  0,  X  *  0  is  analyzed.  The  level  lines  obtained  by  intersecting  this  surface 
by  the  equatorial  plane  of  the  Moon,  the  zero  meridian  plane  (X  =  0),  and  the 
plane  corresponding  the  meridian  X  =  90°  are  presented.  Orig.  art.  has:  4 
figures.  English  translation:  7  pages. 
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DETERMINATION  OF  THE  GRAVITATIONAL  FIELD  OF  THE  MOON 
BY  THE  MOVEMENT  OF  THE  ARTIFICIAL  MOON  SATELLITE  LUNA-10 

E.  L.  Akim 

An  examination  is  made  of  the  movement  of  the  artificial  satellite 
of  the  moon,  Luna-10,  in  the  central  field  of  gravitation  of  the  moon 
with  the  taking  into  account  of  the  gravitational  effect  of  the  earth, 
the  sun.  and  the  planets. 

If  the  gravitational  field  of  the  moon  were  central  and  the 
effect  of  external  bodies  did  not  exist  the  motion  of  the  satellite 
would  be  along  an  undisturbed  orbit  (the  Kepler  ellipse),  the  form  and 
dimensions  of  which  would  remain  unchanged  in  absolute  space.  The 
noncentral  quality  of  the  gravitational  field  of  the  moon  and  the 
action  of  the  outer  bodies,  the  principal  ones  of  which  are  the  earth 
and  the  sun  bring  about  disturbances  in  the  motion  of  the  satellite. 
Under  the  action  of  the  disturbances  of  the  orbit,  the  satellite  in 
the  course  of  time  evolutionizes. 

The  disturbing  influence  of  the  earth  and  the  sun  on  the  movement 
of  the  satellite  of  the  moon  is  well  known.  The  greatest  interest 
is  afforded  by  the  evolution  of  the  orbit  of  the  satellite  in  respect 
of  the  unknown  noncentral  quality  of  the  field  of  gravitation  of  the 
moon.  Knowledge  of  this  evolution  enables  one  to  determine  the 

FTD-HT-67-50  1 


parameters  of  the  gravitational  field  of  the  noon. 

For  revealing  the  evolution  of  the  orbit  of  the  artificial 
satellite,  Luna-10,  recourse  was  had  to  measurements  of  the  trajectory 
of  the  movement  of  the  satellite  occurring  during  the  course  of  the 
whole  time  of  Its  active  existence  (3rd  of  April  to  the  30th  of  May, 
1966).  The  trajectory  measurements  were  submitted  to  statistical 
processing  for  the  purpose  of  joint  determination  of  the  parameters 
of  the  gravitational  field  of  the  moon  and  the  elements  of  the  orbit 
of  the  satellite.  At  the  basis  of  the  procedure  of  the  processing 
of  the  measurements  there  has  be.en  put  the  analytical  theory  of  the 
movement  of  the  satellite  of  the  moon,  which  made  it  possible  to 
embrace  by  a  single  computation  the  whole  two-month  measuring  interval 

of  the  motion  of  the  satellite. 

*  » 

.  The  description  of  the  selenocentric  motion  of  the  satellite  and 
the  motion  of  the  moon  around  Its  own  center  of  mass  was  worked  out 
in  the  Cartesian  rectangular  selenocentric  system  of  coordinates  XYZ. 
The  plane  XY  of  this  system  corresponds  with  the  plane  of  the  average 
equator  of  the  moon,  and  the  plane  XZ  with  the  plane  of  Its  zero 
meridian  of  the  period  t^.  The  directions  of  the  axes  of  the  system 
are  motionless  with  respect  to  the  stars.  The  axis  X  of  the  system 
is  directed  to  the  earth,  the  axis  Z  to  the  North  Pole  of  the  moon, 
and  the  axis  Y  completes  the  system  to  the  right.  For  the  description 
of  the  motion  of  the  satellite  there  are  used  the  elements  of  its 
orbit;  the  great  semiaxis  a,  the  eccentricity  e,  the  inclination  1, 
the  longitude  of  the  rising  node  ft,  the  angular  distance  of  the  peri- 
center  from  the  node  o>,  and  the  time  of  the  passing  of  the  node  Tq. 
(The  reading  off  of  the  angular  orbital  elements  i,  ft,  o>  is  done  from 
th  niane  of  the  average  equator  of  the  moon  and  its  zero  meridian 
of  the  period  tQ  by  the  method  accepted  in  celestial  mechanics).  It 
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Is  presumed  that  the  turning  of  the  moon  around  its  own  center  of 
mass  occurs  in  accordance  with  Cassini’s  laws  —  uniformly  around  a 
motionless  axis  OZ  of  the  introduced  system  of  coordinates.  The 
motion  of  the  center  of  the  mass  of  the  moon  (beginning  of  the  system 
of  coordinates)  in  the  geoequator ial  system  of  coordinates  with 
average  equinox  of  the  period  of  i960  (0)  is  given  by  Brown's  theory. 

The  noncentral  quality  of  the  gravitational  field  of  the  moon  is 
an  essential  fact  determining  the  evolution  of  the  satellite  Luna-10. 

The  disturbances  of  the  orbit  of  the  satellite  which  arise  because  of 
the  noncentral  quality  of  the  field  of  gravitation  of  moon  are  revealed 

with  special  clarity  in  the 
evolution  of  the  longitude  (2 
of  the  rising  node  of  the 
satellite  and  the  angular 
distance  <d  of  the  pericenter 
from  the  node.  This  evolu¬ 
tion  of  the  elements  (2  and  o> 
of  the  orbit  of  the  satellite 
during  the  time  of  ita  active 
exlstance  is  3hown  in  Fig.  1 
as  a  function  of  the  number 
of  turns  of  the  satellite. 

The  evolution  of  the  param¬ 
eters  (2  and  co  contains  an 
obvious  secular  drift  leading  to  a  regression  of  both  parameters.  Onto 
the  secular  drift  of  the  parameter  o>  there  is  superimposed  a  noticeable 
periodical  disturbance.  For  460  turns  of  the  satellite  the  disturbances 
of  the  parameters  (2  and  co,  brought  about  by  the  noncentral  quality  of 
the  field  of  gravitation  of  the  moon  reach  the  values  Afl  «  -7. 7° 


Fig.  1.  Disturbances  in  longitude  ft  of 
the  rising  node  (unit)  and  angular 
distance  of  the  pericenter  from  the  node 
of  the  orbit.  The  heavy  lines  mean  on 
account  of  the  noncentral  quality  of  the 
field  of  gravitation  of  the  moon  and  the 
thin  lines  mean  on  account  of  the  gravi¬ 
tational  influence  of  the  earth  and  the 
sun. 
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and  Au>  =  -11.8°.  The  disturbances  of  the  inclination  i  of  the  orbit 
of  the  satellite  and  its  eccentricity  e  have  basically  a  periodical 
character  and  have  an  amplitude  of  Ai  «  0.15  and  Ae  «  0.005. 

The  disturbances  of  the  orbit  of  the  satellite  on  account  of  the 
noncentral  quality  of  the  field  of  gravitation  of  the  moon  lead  to 
disturbances  of  its  coordinates,  amounting  for  one  turn  of  the  satellite 
to  the  value  |ap}  «  0.75  km. 

The  gravitational  influence  of  the  earth  and  the  sun  on  the 
motion  of  Luna-10  also  leads  to  the  evolution  of  its  orbit.  This 
Influence  brings  about  a  regression  of  the  node  and  peri  center  of  the 
orbit  of  the  satellite,  which  are  also  shown  in  Fig.  1.  During  the 
time  of  the  active  existence  of  the  satellite  these  disturbances  amount 
to  the  value  Aft  =  -1°  and  Atu  =  -2°.  Under  the  influence  of  the  earth 
and  the  sun  the  eccentricity  of  the  orbit  of  the  satellite  becomes 
less  in  the  Interval  of  time  under  consideration.  In  an  almost 
unchangeable  great  semiaxis  of  the  orbit  the  lessening  of  the  eccentric¬ 
ity  leads  to  an  increase  in  the  perlcenter  of  the  orbit.  The 
disturbances  of  the  coordinates  of  the  satellite  Luna-10  brought  about 
by  the  action  of  the  earth  and  the  sun  for  one  turn  of  the  satellite 
do  not  exceed  the  value  [  A? 1  «  0.11  km. 

The  comparative  evaluations  presented  show  that  for  the  orbit  of 
the  Luna-10  the  cisturbances  brought  about  by  the  noncentral  quality 
of  the  field  ux  gravitation  of  the  moon  by  a  factor  of  5  or  6  exceed 
the  disturbances  brought  about  by  the  gravitational  effect  of  the 
earth  and  the  sun.  The  planet  disturbances  of  the  motion  of  the 
satellite  in  the  interval  of  time  under  consideration  are  small  and 
are  therefore  not  taken  into  the  reckoning. 

The  processing  of  the  trajectory  measurements  containing  their 
sinpJe  dynamic  bond  on  the  interval  of  flight  of  the  satellite  Luna-10 
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in  two  lunar  months  made  it  possible  to  determine  the  quantitative 
characteristics  of  the  noncentral  quality  of  the  gravitational  field 
of  the  moon. 

The  expression  for  the  gravitational  potential  U  of  the  moon  is 
taken  in  the  form  of  an  expansion  into  series  in  accordance  with  the 
spherical  functions 

*)--£{* +2  s(-r)"x' 

'  It-*  M 

[cBm  coamX+Oin  mk]P„n(m\ 

where  u  is  the  mass  of  the  moon,  R  is  its  average  radius;  r,  and  X 
are  the  spherical  coordinates  of  the  points:  r  is  the  polar  radius, 
f  is  the  latitude  read  from  the  average  equator  of  the  moon,  X  is  the 
longitude  read  from  the  zero  meridian  of  the  period  tQ,  and  Pnm(sin  $) 
are  the  attached  Legendre  functions.  As  the  parameters  to  be  determined 
of  the  gravitational  field  of  the  moon  there  are  used  the  coefficients 
cnm  and  dnm  of  this  analysis. 

As  the  result  of  the  processing  of  the  trajectory  measurements, 
there  are  obtained  numerical  values  for  eleven  coefficients  of  the 
analysis  of  the  gravitational  potential  of  the  moon.  These  values 
together  with  'heir  maximum  possible  errors  are  given  below. 

c20  "  (-°-2°6  2  0.022)  •  10-5, 
c21  -  (0.157  +  0.059)  *  10"4, 
d21  =  (0.561  +  0.558)  ’  10-5, 
c22  =  (0.140  +  0.012)  *  10“4, 
d22  -  (-0.159  +  0.145)  •  10"5, 

c50  “  (-0. 565  +  0.099)  •  10-4, 

c5:L  -  (-0.568  +  0.026)  •  10“4, 

d51  »  (-0.178  +  0.052)  *  10-^, 
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Fig.  2.  Sections  of  the 
surface  of  tne  level  of 
gravitational  potential  of 
the  moon  (radial  inclina¬ 
tion  from  the  circumference 
increased  by  1,000  times). 

A  is  the  equatorial,  £  the 
meridional  (A  =  0),  and  B 
is  the  meridional  (A  =  90°). 
KEY:  1)  In  the  direction 
of  the  earth. 


c-22  =  (0.118  +  0.047)  •  10"^, 
d32  =  (-0.702  +  4.595)  •  10-6, 
c4o  =  (0.25}  +  0.270)  *  iO-1*. 

Between  the  errors  in  the  determina¬ 
tion  of  the  parameters  C2o  and  C4q  there 
is  a  substantial  correlation.  The 
coefficient  of  the  correlation  k  =  O.99. 
The  mutual  correlation  between  the  errors 
of  the  remaining  parameters  to  be 
determined  does  not  go  above  the  figure 
0.4. 

The  numerical  values  presented  for 
the  parameters  C20  and  C22  aSree  with 
the  values  known  for  them  obtained  in 
accordance  with  the  libration  measure¬ 
ments  (1). 

The  difference  from  zero  of  the 
coefficients  cnm  with  the  odd  suffix  m 
and  the  coefficients  dnm  with  the  even 
suffix  m  is  evidence  that  the 
gravitational  fields  on  the  side  of  the 
moon  visible  from  the  earth  and  on  that 
invisible  from  the  earth  are  not 
symmetrical. 

For  illustration  of  the  determined 
gravitational  potential  of  the  moon, 
there  is  considered  the  surface  of  a 
plane  running  through  a  point  with  the 
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spherical  coordinates  r  =  1738  km,  ^  =  0,  X  =  0.  In  the  Pig.  2  there 
are  constructed  lines  of  the  plane  or  level  obtained  by  intersections 
of  this  surface  by  the  equatorial  plane  XY  (Pig.  2  A),  the  plane  of  the 
zero  meridian  XZ  (Fig.  2  6),  and  the  meridional  plane  YZ  corresponding 
to  the  longitude  X  =*  90°  (Pig.  2  B).  The  radius  of  the  circumference 
in  these  drawings  corresponds  to  the  average  radius  of  the  moon,  equal 
to  1738  km.  The  radial  deviation  of  the  lines  of  the  level  from  the 
circumference  for  clarity  is  drawn  with  an  exaggeration  of  a  thousand 
times. 

There  is  seen  the  pear-shaped  quality  of  the  surface  of  the 
potential  with  the  bulging  on  the  reverse  side  of  the  moon.  The  pear- 
shaped  effect  is  formed  basically  under  the  influence  of  terms  in  the 

l 

analysis  of  the  potential  which  contain  the  tesseral  harmonic  P^ 

(with  the  coefficient  c^)  and  the  zonal  harmonic  P2°. 

The  results  presented  above  are  only  preliminary.  Further 
processing  of  the  trajectory  measurements  in  accordance  with  Luna-10 
and  analysis  of  the  movement  of  succeeding  satellites  of  the  moon  will 
make  it  possible  to  refine  and  supplement  the  parameters  obtained  for 
the  gravitational  potential  of  the  moon  and  also  to  refine  the  movement 
of  its  center  of  mass. 

Entered 

September  13,  1956 
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